Extension of the spin-1/2 frustrated square lattice model: 
the case of layered vanadium phosphates 
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We study the influence of the spin lattice distortion on the properties of frustrated magnetic 
systems and consider the applicability of the spin-1/2 frustrated square lattice model to materials 
lacking tetragonal symmetry. We focus on the case of layered vanadium phosphates AA'VO(P04)2 
(AA' = Pt)2, SrZn, BaZn, and BaCd). To provide a proper microscopic description of these com- 
pounds, we use extensive band structure calculations for real materials and model structures and 
supplement this analysis with simulations of thermodynamic properties, thus facilitating a direct 
comparison with the experimental data. Due to the reduced symmetry, the realistic spin model 
of layered vanadium phosphates AA'VO(P04)2 includes four inequivalent exchange couplings: Ji 
and J[ between nearest-neighbors and J2 and J2 between next-nearest-neighbors. The estimates 
of individual exchange couplings suggest different regimes, from J[/Ji and J2/J2 close to 1 in 
BaCdVO(P04)2, a nearly regular frustrated square lattice, to J[/Ji — 0.7 and J2/J2 — 0.4 in 
SrZnVO(PC>4)2, a frustrated square lattice with sizable distortion. The underlying structural differ- 
ences are analyzed, and the key factors causing the distortion of the spin lattice in layered vanadium 
compounds are discussed. We propose possible routes for finding new frustrated square lattice mate- 
rials among complex vanadium oxides. Full diagonalization simulations of thermodynamic properties 
indicate the similarity of the extended model to the regular one with averaged couplings. In case of 
moderate frustration and moderate distortion, valid for all the AA'VO(P04)2 compounds reported 
so far, the distorted spin lattice can be considered as a regular square lattice with the couplings 
(Ji + Ji)/2 between nearest-neighbors and (J2 + J2)/2 between next-nearest-neighbors. 



PACS numbers: 75.10.Jm, 75.30.Et, 75.50.-y, 71.20.Ps 
I. INTRODUCTION 

Frustrated spin systems represent one of the actively 
developing topics in solid state physics. The vast in- 
terest in magnetic frustration originates from a num- 
ber of unusual phenomena (spin-liquid ground stated the 
formation of supersolid phases in high magnetic fields, 2 
etc.) suggested by theory. The theoretical predictions 
challenge the experiment that, however, requires proper 
frustrated materials. The search for the respective com- 
pounds has been a long story in inorganic chemistryP^ 
and the problem turned out to be quite complex. To 
meet the theoretical predictions, one has to find a mate- 
rial that reveals frustrated geometry of magnetic atoms 
and presents spin degrees of freedom only to avoid any 
foreign effects (e.g., orbital ordering) tending to lift the 
frustration. A number of frustrated spin models still 
lack the proper realizations, especially for the case of 
spin-1/2, where the strongest quantum effects and the 
most interesting phenomena are expected. For other 
models, few appropriate materials are known and exten- 
sively studied. For example, the mineral herbertsmithite 
ZnCu3(OH) 6 Cl2 was recently proposed as a spin-1/2 
kagome material P However, the actual physics of this 
compound is still debated due to the Cu/Zn antisite dis- 
order and the presence of non-magnetic sites within the 
kagome layersP Other natural kagome materials - min- 
erals kapellasite and haydeeite - are proposed. Their 
structures do not suffer from the disorder effects, but the 
pure kagome physics is again modified due to the non- 



negligible interactions beyond nearest neighbors.^ 

Clearly, the task of finding an ideal frustrated material 
is hardly solvable at all. Therefore, it is instructive to ex- 
amine which deviations from the ideal spin model may 
be allowed and, to a certain extent, do not qualitatively 
modify the properties of this model. Considering struc- 
tural distortions and the resulting spin lattice distortions 
is especially attractive, since lots of known materials have 
low symmetry in contrast to the high geometrical sym- 
metry of the most theoretically studied frustrated spin 
models. 

The major part of the frustrated magnetic materials 
arc the so-called geometrically frustrated magnets. In 
these systems, the frustration arises due to the competi- 
tion of equivalent exchange couplings: in the most sim- 
ple case, antiferromagnetic (AFM) couplings on a trian- 
gle. Then, any structural distortion should inevitably 
reduce this competition hence reducing the frustration. 
In the other group of the frustrated materials, the com- 
peting interactions are inequivalent but their topology 
and magnitudes can be tuned so that the strong quan- 
tum fluctuations destroy the long-range ordering, simi- 
lar to the geometrically frustrated magnets. The latter 
group looks more favorable to tolerate the structural dis- 
tortions, since the modification of the spin lattice can 
probably be balanced by the ratios of the competing in- 
teractions hence preserving the strong frustration. To 
study this issue in more detail, we focus on a specific 
model - spin-1/2 frustrated square lattice (FSL) - and 
consider a number of recently proposed FSL materials. 
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FIG. 1: (Color online) Phase diagram of the FSL modePand 
the respective model compounds (see text for references). The 
inset shows the regular FSL with the NN (Ji) and NNN (J 2 ) 
couplings denoted by solid and dashed lines, respectively. 



The ideal (regular) FSL model assumes two competing 
interactions: the nearest-neighbor (NN) interaction J\ 
running along the side of the square and the next-nearest- 
neighbor (NNN) interaction J2 running along the diago- 
nal of the square (see the inset of Fig. [T]) . The phase dia- 
gram (Fig. [lj reveals three ordered phases (ferromagnet, 
Neel antiferromagnet, and columnar antiferromagnet) 
and two critical regions around J2/J1 — ±0.5, where a 
spin-liquid ground state is expected Recent theoretical 
studies also considered the extended model with inequiv- 
alent NN couplings J x and J[ (Refs. ll0lllll2ll3ll4l) . This 
spatial anisotropy tends to narrow the critical region^ 
and to destroy it completely at a certain value of J[/ J\. 
Therefore, the distortion is not favorable for the frustra- 
tion, but the specific geometry of the lattice enables to 
preserve the strong frustration and the resulting critical 
region at moderate distortion. Below (Sec. VI), we will 
show that the anisotropy of the NNN couplings on the 
square lattice should have an even weaker (and, likely, 
opposite) effect on the frustration. 

Experimental studies of the FSL systems have utilized 
a number of model compounds. First reports focused 
on Li2VOXC>4 (X = Si, Ge) materials that lay far away 
from the critical regions and did n ot show any sp ecific 
properties caused by the frustration! 16 * 17 * 18 * 19 * 20 * 21 ! Later 
on, the vanadyl molybdate VOM0O4 was found to reveal 
unusual structural changes upon cooling, and this effect 
was tentatively ascribed to the magnetic frustration de- 
spite the fru stratio n ratio J2/J1 was quite small (likely, 
J2/J1 < 0.2).l22122lThe region of ferromagnetic (FM) J x 
- AFM J2 was accessed by studying the layered vana- 
dium phosphate Pb 2 VO(P0 4 ) 2 (Refs. 124125126]) and the 
relate d AA / VO(P0 4 )2 compounds with AA' = SrZn and 
BaZn. 25 27 Quite recently, we proposed two more com- 
pounds, PbV0 3 (Refs. EE] and dj) and BaCdVO(P0 4 ) 2 



(Ref. [30J) 7 that lay very close to the critical regions at 
J2/J1 = 0.5 and —0.5, respectively. Interestingly, the 
latter material lacks the tetragonal symmetry and reveals 
a distorted FSL. Nevertheless, we succeeded to observe 
two effects predicted for the regular FSL: the suppres- 
sion of the specific heat maximurrP and the pronounced 
bending of the magnetization curved 

In other systems, the problem of the spin lattice dis- 
tortion may also be crucial. The layered copper oxy- 
chloride (CuCl)LaNb2C>7 was recently proposed as a 
promising FSL material, lacking long-range magnetic 
orc j er |32j u owever careful studies indicated the struc- 
tural distortion^! that completely changed the under- 
lying spin model and precluded from any interpreta- 
tions within the FSL frameworkP3 Thus, it is impor- 
tant to understand whether the FSL model can be ap- 
plied to BaCdVO(PC>4)2 and, more generally, to the low- 
symmetry materials. Below, we study this issue in de- 
tail, discuss the whole family of the AA'VO(P04)2 com- 
pounds (hereinafter, we imply that AA' = Pb2, SrZn, 
BaZn, and BaCd), and provide quantitative estimates 
for the spin lattice distortion. Our approach combines 
several computational methods (band structure calcula- 
tions, subsequent analysis of the exchange couplings, and 
model simulations) in order to analyze magnetic interac- 
tions in these systems, derive the proper spin model, and 
facilitate the comparison with the experimental data. We 
show that the structural distortion in the AA'VO(P04)2 
compounds is a minor effect as compared to the frustra- 
tion, and the FSL description holds. 

The outline of the paper is as follows. We start with an 
analysis of the crystal structures in Sec. Inland review the 



computational methods employed in our work (Sec. Ill I. 
Then, we address several problems: (i) the realistic spin 
model and the magnitude of the distortion (Sec. IV I; (ii) 
structural factors that influence the distortion of the spin 
lattice (Sec. Vj); and (iii) ther modyna mic p roperties of 



VII we discuss 



the extended model (Sec. |VI[ ). In Sec. 

the results of our study that suggests an accurate (and, 
within the Heisenberg model, exact) way to treat the 
FSL-like spin systems of the AA'VO(P04)2 phosphates. 
We also present a more general recipe for finding strongly 
frustrated square lattices in the compounds with similar 
topology of the magnetic layer, being generic for most of 
the FSL materials. Finally, we show how the distortion 
of the square lattice affects thermodynamic properties of 
the model and the magnitude of the frustration. 



II. CRYSTAL STRUCTURES AND 
EXPERIMENTAL RESULTS 

Most of the vanadium-based FSL compounds reported 
so far reveal similar [VOXO4] magnetic layers. These 
layers are built of VO5 pyramids and XO4 tetrahedra 
(see left panel of Fig. [2]) with X being a non-magnetic 
cation (P, Si, Ge, or Mo +6 ). The connections via the 
tetrahedra provide superexchange pathways for both NN 
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FIG. 2: (Color online) Crystal structure of the AA'VO(P04)2 compounds and the underlying spin model. The left panel shows 
a single [VOPO4] layer. The middle panel presents the stacking of the [VOPO4] layers and the [AA'P04] blocks as well as 
the angle ip measuring the layer buckling: larger and smaller spheres denote the A and A' cations, respectively. The right 
panel shows the magnetic interactions in the [VOPO4] layers (compare with the regular model in the inset of Fig. [I]): solid, 
dash-dotted, dashed, and dotted lines indicate Ji, J[, J 2 , and J 2 , respectively. 



and NNN couplings. Vanadium atoms have the steady 
oxidation state of +4 (implying the electronic configura- 
tion 3d 1 and spin-1/2), while the valence of the X cation 
can be changed and controls the filling of the interlayer 
space. Thus, the interlayer space is empty in VOM0O4 
(Ref.EU), filled by lithium atoms in Li 2 VOX0 4 (X = Si, 
Ge^PH and filled by complex inter layer [A A'PC^] blocks 
in the AA'VO(P0 4 ) 2 phosphates! 37 ' 38 ' 391 The magnetic 
layer is compatible with the tetragonal symmetry, hence 
the FSL model can be realized. Yet the overall symmetry 
of the structure is sometimes reduced due to the complex 
configuration of the interlayer block, and this is the case 
for the AA'VO(P0 4 )2 compounds. 

The crystal structure of the AA'VO(PC>4)2 phosphates 
(Fig. |2| is fairly flexible with respect to the accommo- 
dation of different divalent metal cations A and A' in 
the interlayer block. Most of the compounds combine a 
larger (Ba or Sr) and a smaller (Cd or Zn) cation, yield- 
ing the ortho rhom bic symmetry of the structure (space 
group Pbca) ! 37 ' 38 l Yet it is also possible to use the same 
cation (Pb) for A and A', and the resulting structure 
has monoclinic symmetry (space group c)P Differ- 
ent symmetries produce slightly different geometries of 
the magnetic [VOPO4] layers, hence different types of 
the spin lattice distortion should be expected. In the fol- 
lowing, we will focus on the magnetic layer typical for 
the orthorhombic AA'VO(PC>4)2 compounds with dif- 
ferent A and A' (Fig. [21. To ease the comparison to 
monoclinic Pb2VO(PC>4j2, one should convert the crys- 
tal structure^ to the non-standard setting P2i/b with 
a being the monoclinic axis (in this setting, the layer 
corresponds to one shown in Fig. 111. Basically, the 
whole discussion of the spin lattice distortion is applica- 
ble to Pb2VO(PC>4)2, but due to the lower lattice sym- 
metry this compound has one additional feature: two 
incquivalcnt NN couplings along the b axis. However, 
our estimates show that these couplings nearly match. 
Therefore, the spin system of Pb2VO(PC>4)2 can be ad- 
equately described with four parameters similar to the 



other AA'VO(P0 4 ) 2 compounds. 

The lack of the tetragonal symmetry in the 
AA'VO(PC>4)2 phosphates gives rise to four different in- 
teractions in the magnetic [VOPO4] layers. According 
to Fig. [2] we label the NN interactions as Ji, J[ and 
the NNN interactions as Ji,J'i- Previous experimen- 
tal works implicitly assumed the ideal FSL model with 
J\ = J[ and J2 = J2 as a na tural, albeit therein un- 
justified, approximation. ^ ' 25 ' 30 ' In the previous studies, 
experimental data on the magnetic susceptibility and the 
specific heat were fitted with high-temperature series ex- 
pansions (HTSE) for the regular FSLA^to yield the effec- 
tive couplings Ji xp and J| xp (Table[j). The FM J x - AFM 
Ji regime of the FSL was further supported by analyzing 
field dependence of the magnetization of BaCdVO(P04)2 
(Ref. |30|) and the ground states of Pb 2 VO(P04)2 and 
SrZnVO(P0 4 ) 2 ( Ref s. [Ml and [57]). 

Within the phenomenological approach, the consistent 
interpretation of the experimental results justifies a pos- 
teriori the application of the ideal FSL model to the 
AA'VO(PC>4)2 compounds. Of course, a microscopic ap- 
proach requires a more careful consideration of all the 
four incquivalcnt exchange couplings in the magnetic 
[VOPO4] layers. However, it is quite difficult (at least, 
experimentally) to go beyond the regular FSL descrip- 



TABLE I: Experimental exchange couplings ( Jl xp , J 2 xp ) eval- 
uated within the regular FSL model and the resulting frus- 
tration ratios in the AA'VO(P04)2 compounds 
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tion due to the lack of theoretical results for the extended 
J i — J[ — J 2 — J' 2 model. In the following, we address the 
problem using computational methods. These methods 
are known to provide a reliable microscopic description 
of complex spin systems, including the FSL compounds 
Li 2 VOX0 4 (Refs. HE and [HJ). Yet we also refer to the 
phenomenological results and show that at sufficiently 
high temperatures the distorted FSL can be considered 
as a regular FSL with effective, averaged NN and NNN 
couplings. 



III. METHODS AND MODELING 

Scalar-relativistic band structure calculations were 
performed within the full-potential local-orbital scheme 
(FPLO7.00-27)p^and the exchange-correlation potential 
by Perdew and WangES was applied. Our calculations 
employed experimental crystal structures and a number 
of modified and model structures, as described below. 
Different fc-meshes were used depending on the size and 
the geometry the unit cell. In all the calculations, the 
convergence with respect to the fc-mesh was carefully 
checked. 

To evaluate the exchange couplings in the 
AA'VO(P04)2 compounds, we use two different 
approaches. First, local density approximation (LDA) 
calculations are performed. These calculations enable 
to select relevant states and to estimate hopping pa- 
rameters (t) for the respective bands by fitting these 
bands with a tight-binding (TB) model. The LDA 
calculations fail to reproduce the strong correlation 
effects in the vanadium 3d shell; therefore, the corre- 
lations are included on a model level. The hoppings 
are introduced to an extended Hubbard model with 
the effective on-site Coulomb repulsion U c s — 4.5 eV 
(this value is representative for vanadium oxides, see 
Refs. OS [19J S2J and |43j. The strongly correlated limit 
t <c U c g and the half-filling regime justify the reduction 
to the Heisenberg model for the low-lying excitations. 
Then, AFM contributions to the exchange couplings are 
estimated as J 4 AFM = At^/U c s- Within this approach, all 
the possible AFM couplings are evaluated. 

Second, we consider the correlation effects within the 
self-consistent calculations and employ the local spin den- 
sity approximation (LSDA)+[7 method. Total energies 
for a number of ordered spin configurations are mapped 
onto the classical Heisenberg model to yield the estimates 
for both FM and AFM couplings, hence supplementing 
the TB analysis. LSDA+t/ treats correlation effects in 
a mean-field approximation and uses two input parame- 
ters, Ud and Jd, to describe the on-site Coulomb repulsion 
and the intraatomic exchange, respectively. Since the ex- 
change parameter has minor influence on the results, we 
fix Jd — 1 eV as a representative value. Yet the choice of 
the repulsion parameter Ud may be crucial, especially in 
case the exchange couplings are weak (see Ref. 2Hfor an 
instructive example). To reduce the ambiguity related 



to the choice of Ud, we thoroughly compare the compu- 
tational results with the experimental data, as further 
discussed in Sec. IIVI 

The LDA calculations employed the full symmetry of 
the crystal structures, 112-atom orthorhombic unit cells, 
and the fc-mesh of 256 fc-points with 75 points in the 
irreducible part of the first Brillouinc zone (IBZ). In case 
of Pb2VO(P04)2, the 56-atom monoclinic unit cell and 
a mesh of 512 fc-points (170 in IBZ) were used. 

To realize different spin orderings within the LSDA+C/ 
calculations, one has to reduce the crystal symmetry and, 
in some cases, to extend the unit cell. The evaluation of 
the four exchange couplings in the AA'VO(P04)2 com- 
pounds requires the doubling of the unit cell in the a di- 
rection, hence 224-atom unit cells should be constructed. 
For such unit cells, full-potential calculations are ex- 
tremely time-consuming and, likely, not accurate enough 
for a reliable evaluation of the rather small exchange 
constants.^ Therefore, in our LSDA+C/ calculations 
we simplify the crystal structures of the AA'VO(P04)2 
compounds and construct a number of modified struc- 
tures. The idea resembles our study of Ag2VOP207 
(Ref. l42|l : the leading magnetic interactions take place 
in the V-P-0 layers, hence it is essential to use the cor- 
rect geometry of the layer, while the stacking of the layers 
and the filling of the interlayer space have minor effect 
on the leading exchange couplings. 

To build the simplified structures of the AA'VO(P0 4 ) 2 
compounds, we keep the exact geometry of the [VOPO4] 
layers, stack these layers one onto another, and fill the 
interlayer space with lithium atoms, providing the proper 
charge balance. The resulting composition is LiVOP04. 
The interlayer separation is fixed at 6.5 A to achieve real- 
istic, sufficiently weak interlayer hoppings (below 2 meV). 
To justify the structure simplification, we perform the 
TB analysis. The difference between the respective hop- 
pings in the experimental and modified structures does 
not exceed 5 %, implying an error below 10 % for the 
exchange couplings. Such an error is definitely accept- 
able for the further LSDA+C/ calculations. Basically, 
the structure simplification provides a promising com- 
putational approach to the magnetic properties of low- 
dimensional spin systems with complex crystal struc- 
tures. In our LSDA+f/ calculations, we use 64-atom su- 
percells [2a x b x c{— 6.5 A)] with triclinic symmetry 
(space group PI) and a mesh of 108 fc-points. 

Finally, we also performed a number of LDA calcula- 
tions for model structures in order to study the influence 
of individual structural changes on the spin lattice dis- 
tortion. The model structures were built similar to the 
simplified structures described above. The initial geome- 
try was taken from the structure of ai-LiVOP04 that in- 
cludes regular [VOPO4] layers separated by Li cations.^ 
Then, we introduced a number of structural distortions 
and checked the changes of the exchange couplings (see 
Sec. [V]for details). The initial structure reveals a rather 
low interlayer separation of 4.45 A and yields sizable in- 
terlayer interactions. To reduce the interlayer interac- 
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tions and to properly emulate the two-dimensional (2D) 
character of the FSL compounds, we increased the inter- 
layer spacing up to 6.5 A. We also reduced the symmetry 
down to the orthorhombic space group Pbma [Pbcm in 
the standard setting) that allowed the distortions of the 
[VOPO4] layer. A fc-mesh of 4096 points (729 in the IBZ) 
was used. 

Full diagonalization (FD) simulations were performed 
for the N = 16 (4 x 4) cluster using the ALPS simu- 
lation package.^ Basic thermodynamic quantities (mag- 
netic susceptibility and specific heat) were evaluated by 
an internal procedure of the program. In general, the FD 
simulations suffer from finite-size effects, because current 
computational facilities do not allow to perform the cal- 
culations for large clusters. The presently available clus- 
ter size (normally, 16 or 20 sites) is sufficient to obtain 
the quantitatively correct information on thermodynamic 
properties of one-dimensional spin systems only. The ac- 
curacy of the FD simulations for two-dimensional systems 
is challenged by the experi mental da ta and the results of 
other simulation techniques! 25 * 30 * 4 ^ Yet, the FD simula- 
tions are able to provide qualitatively correct trends upon 
the change of the model parameters (e.g., the change of 
the frustration ratio in the FSL, see Ref. [9j). Keeping 
in mind these considerations, we restrict ourselves to the 
analysis of the relative changes in the thermodynamic 
properties, as the spin lattice is distorted. Investigation 
of the ground state and quantitative simulation of ther- 
modynamic properties for the extended FSL model are 
clearly beyond the scope of the present paper that in- 
tends to stimulate further studies of the problem. 



IV. EXCHANGE COUPLINGS IN AA'VO(P0 4 ) 2 
A. LDA and tight-binding analysis 

In this section, we consider the band structures of the 
AA'VO(P04)2 compounds and evaluate individual ex- 
change couplings. The LDA density of states (DOS) for 
BaCdVO(P0 4 ) 2 is shown in Fig. [3] This plot is rep- 
resentative for the whole family of the AA'VO(P04)2 
materials. The electronic structure resembles that of 
other vanadium phosphates P^l The states below —2.5 eV 
are mainly formed by oxygen orbitals, while the states 
close to the Fermi level reveal predominant vanadium 
contribution. The vanadium bands are rather narrow 
and show 3d-related crystal field levels (see the inset of 
Fig. [3| as expected for the square- pyramidal or distorted 
octahedral coordination of V +4 (Ref. |4"9"]| . The lowest- 
lying vanadium states formed by d xy orbitals lie at the 
Fermi level. The respective orbitals are located in the 
basal planes of the square pyramids hence overlapping 
with p orbitals of the basal oxygen atoms and facilitat- 
ing exchange couplings in the ab plane. The LDA energy 
spectrum is metallic in contradiction to the experimen- 
tal green or yellow-green color that indicates insulating 
behavior. The lack of the energy gap is a typical fail- 
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FIG. 3: (Color online) LDA density of states for 
BaCdVO(P04)2 (the contribution of barium is not shown, 
because it is negligible in the whole energy range) . The Fermi 
level is at zero energy. The inset shows the orbital resolved 
DOS for vanadium. 



ure of LDA due to the underestimate of strong electron 
correlations in the 3d shell. Applying LSDA+J7, we find 
insulating spectra with an energy gap of 1.5 — 2.8 eV 
depending on the Ud value. The upper estimate is in 
reasonable agreement with the experimental sample col- 
ors. 

For the TB analysis, we select the eight half-filled d xy 
bands lying at the Fermi level (Fig. |4|. These bands 
originate from eight vanadium atoms in the unit cell of 
AA'VO(P0 4 ) 2 . The unit cell includes two [VOPO4] lay- 
ers (see middle panel of Fig. [2]) and four vanadium atoms 
from each of the two layers. The layers are well separated 
by the [AA'P04] block; therefore, interlayer interactions 
are very weak, and the bands are close to double degen- 
eracy in the whole Brillouin zone. The leading interlayer 
hopping is about 1 mcV, implying the AFM interactions 
of about 0.01 K, well below the in-layer interactions (see 
Tables|l]and|n|). In Fig.g we show the band structures of 
BaCdVO(P07) 2 and SrZnVO(P0 4 ) 2 at the Fermi level. 
The two plots are rather similar in the overall behavior 
and in the energy scale, although notable differences are 
found near the X, Y, and T points. The TB analysis 
suggests that the differences are mainly related to the 
NNN couplings £ 2 and t' 2 . 

In Table [IT] we list the leading hopping parameters for 
all the AA'VO(P0 4 ) 2 compounds along with the result- 
ing J AFM values for the two representative and extreme 
cases, SrZnVO(P0 4 ) 2 and BaCdVO(P0 4 ) 2 . All the hop- 
pings beyond t 2 and t' 2 are negligible (below 2 meV) . We 
find a notable distortion of the square lattice in all the 
materials under study, yet the magnitude of the distor- 
tion is rather different. The NNN couplings are AFM 
(see Table p} and correspond to long V-V separations, 
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FIG. 4: (Color online) LDA band structure (thin light lines) and the fit of the tight-binding model (thick green lines) for 
BaCdVO(P04)2 (left panel) and SrZnVO(P04)2 (right panel). The Fermi level is at zero energy. The notation of the k points 
is as follows: F(0,0,0), X (0.5, 0,0), M(0.5, 0.5, 0), Y(0, 0.5,0), Z(0, 0,0.5), and T(0.5, 0, 0.5) (the coordinates are given along 
k x , k y , and k z in units of the respective reciprocal lattice parameters) . 



hence the FM contributions are expected to be small, 
and one could consider the J 2 AFM / J AFM ratio as a good 
estimate for the distortion ratio, J 2 /J 2 , with respect 
to the ideal value J' 2 / J2 — 1 for the regular FSL. We 



find the least pronounced distortion (J 2 \ afm /J£ lFM 
0.83) in BaCdVO(P0 4 ) 2 and the strongest distortion 
(J 2 AFM /J 2 AFM = 0.37) in SrZnVO(P0 4 ) 2 . 

The NN couplings are ferromagnetic (see Table [I]), 
hence the TB model does not yield the direct estimate 
of J'il J\. However, one can assume similar FM contribu- 
tions to J 1 and J[ due to similar V-V separations. Then, 
the difference between the NN couplings should origi- 
nate from the difference between J AFM and J[ AFM , and 
the value of J| AFM — J AFM is the most convenient char- 
acteristic of the distortion. Both J AFM and Jj AFM are 
well below 1 K in BaCdVO(P0 4 ) 2 . In SrZnVO(P0 4 ) 2 , 
J( AFM is about 3 K, while J AFM is still below 1 K. Thus, 
one can expect the sizable difference between J\ and J[ 



TABLE II: The hopping parameters ti (in meV) and the re- 
sulting magnitude of the distortion (J^ A ™/J^ FM ) for all the 
AA'VO(P04)2 compounds along with the antiferromagnetic 
contributions to the exchange integrals J, AFM (in K) for AA' 
= BaCd and SrZn 
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Pb 2 


4" 


15 


47 


38 


0.67 


BaZn 


-11 


12 


46 


36 


0.61 


SrZn 


4 


17 


43 


26 


0.37 




jAPM 


j/AFM 


rAFM 
^2 


j/AFM 


j/AFM jAFM 


BaCd 


0.05 


0.3 


21.0 


17 A 


0.25 


SrZn 


0.2 


3.0 


19.0 


7.0 


2.8 



in SrZnVO(P0 4 ) 2 . In summary, we find weak distor- 
tion of the square lattice in BaCdVO(P0 4 ) 2 and a more 
pronounced distortion for AA' = BaZn and Pb 2 . The 
spin lattice of SrZnVO(P0 4 ) 2 is distorted with respect 
to both the NN and NNN couplings. Below, we will con- 
firm this conclusion with the LSDA+J7 calculations (see 
TablellTTl. 



The TB estimates are in reasonable agreement with 
the experimental data. We find the small NN hoppings 
consistent with the FM nature of J^ xp . On the other 
hand, the sizable hoppings t 2 and t' 2 imply AFM NNN 
couplings in agreement with the AFM coupling J 2 xp . We 
further confirm the FM Ji,J[- AFM J 2 , J 2 scenario with 



"The averaged NN coupling along the b axis, see Sec. [Iljand Ref.l50l 



the LSDA+iy calculations (see Table IIP that perfectly 
match the results of the TB analysis. The remarkable 
agreement of the TB and LSDA+[7 results ensures suffi- 
cient accuracy of our approach. However, the calculated 
values of J 2 and J 2 are likely overestimated as compared 
to J 2 xp , the averaged NNN coupling (Table |l|. We can 
speculate about several reasons for this discrepancy. 

First, one should not expect very precise results of 
band structure calculations while analyzing weak ex- 
change couplings in AA'VO(P0 4 ) 2 . The calculated 
J values likely include a systematic error that, how- 
ever, does not invalidate any of the qualitative conclu- 
sions: not e the studies of the Li 2 VOX0 4 (X = Si, Ge) 
compound d 18 * 19 ! as an instructive example of the correct 
microscopic scenario, emerging from LDA calculations. 
Second, exchange couplings are highly sensitive to the 
geometry of superexchange pathways (see, e.g., Ref. H2J), 
hence it is essential to use the accurate structural in- 
formation in order to obtain precise computational esti- 
mates. In case of AA'VO(P0 4 ) 2 , the crystal structures 
are solved from single-crystal diffraction data, and the 
structural information should be precise. Yet, there are 
some unresolved issues (e.g., poorly reproducible super- 
structure reflections for Pb 2 VO(P0 4 ) 2 , see Ref. l39]l that 
can be crucial for the accuracy of the computational esti- 
mates. Along this paper, we mainly focus on the qualita- 
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17, (eV) I/,(eV) I/„(eV) 

FIG. 5: (Color online) Exchange couplings in Li2VOSi04 (left panel) and Li2VOGeC>4 (middle panel) calculated for different 
values of Coulomb repulsion parameter Ud and the resulting frustration ratios J2/J1 for both the compounds (right panel). 
Shaded stripes show experimental estimates from Refs. .18 19: 2U1 and 1251 



tive differences between the AA'VO(P04)2 compounds, 
while the quantitative analysis should likely reference to 
the experimental data (see Sec. VII for further discus- 
sion) and may require additional structural studies. 



B. LSDA+(7 results 

To get a direct estimate of FM NN couplings in 
AA'VO(P0 4 ) 2 , we turn to the LSDA+J7 approach. In 
this approach, one has to select an appropriate value 
of U d, t he on-site Coulomb repulsion parameter (see 



Sec. trip A number of previous works have estab- 
lished Ud = 3.5 — 4 eV for several V +4 -containing 
compounds j 52 l 53 l S H although the higher value of Ud = 
6 eV^as well as a lower value of U d = 2.3 eV^l were pro- 
posed for specific materials. The reason for the higher Ud 
value in Ref . |55] was tentatively ascribed to the reduced 
p—d hybridization. Still, there is one more uncertainty in 
the selection of this number. The constrained LDA pro- 
cedure and the further LSDA+C/ calculations are usually 
performed within the LMTO-ASA (linearized muffin-tin 
orbitals, atomic spheres approximation) approach that 
employs muffin-tin orbitals. In our work, we use a differ- 
ent basis set (atomic-like local orbitals of FPLO) , hence 
the Ud value may also be different, because the Ud poten- 
tial is applied to the orbitals from the basis set and de- 
pends on the particular choice of the d functions. Indeed, 
our previous FPLO studies showed that the Ud of 5 or 
even 6 eV was required to reproduce the magnetic inter- 
actions in a number of V +4 -containing phosphates! 42 * 44 * 
Keeping in mind the ambiguity of the Ud choice due to 
the differences in the computational method (basis set) 
and in the structural features (different p — d hybridiza- 
tion), we do not use any of the previously established 
Ud values. Rather, we perform calculations for a broad 
range of Ud values and use well-studied and structurally 
similar vanadium compounds as a reference. 

According to Sec. [TTJ the structure of the magnetic 
layer in the AA'VO(P04)2 phosphates is very similar to 
that of Li 2 VOX0 4 compounds (X = Si, Ge). Exchange 
couplings in the latter materials are firmly established by 



fitting ma gnetic susc eptibility and specific heat data with 
the HTSEPI±2120|25| We calculate the exchange couplings 
in Li2VOX04 for a wide range of Ud values and compare 
the results with the experiment (Fig. 5j. In the left and 
middle panels of Fig. [5j one can see that any reasonable 
value of Ud yields the correct energy scale, but there is 
no unique Ud value yielding accurate estimates of both 
J 1 and J 2 . At low Ud, J 2 is overestimated, while high 
Ud values tend to overestimate J\. As long as we are 
interested in the frustration, the essential quantity is the 
frustration ratio, J2/ 1 J\. This quantity is found with a 
sizable error bar due to the uncertainty for the low Ji 
values. Computational estimates of J2/J1 are also quite 
uncertain and highly sensitive to the Ud value (see right 
panel of Fig. [5j. At Ud > 3 eV, we find very low frustra- 
tion ratios, contradicting the J 2 ^S> J\ regime established 
experimentally. The narrow range of Ud = 2 — 3 eV is 
able to reproduce the reasonable frustration ratios, while 
the Ud values below 2 eV lead to FM J\. Thus, we argue 
that the Ud value of 2 — 3 eV should be optimal for repro- 
ducing exchange couplings in the layered vanadium FSL 
compounds. Yet one can construct the correct picture 
solely based on the LSDA+J7 results only in a narrow 
range of the Ud values. 

In the LSDA+f/ study, we restrict ourselves to 
BaCdVO(P0 4 ) 2 and SrZnVO(P0 4 ) 2 as the "edge" mem- 
bers of the AA'VO(P04) 2 series: these compounds show 
the least and the most pronounced distortion, respec- 
tively (see Table |n|. The LSDA+C/ results are summa- 
rized in Table |III| The exchange couplings show a sizable 
dependence on Ud, even in the narrow range of Ud = 2 — 3 
eV. Nevertheless, we find the experimentally observed 
FM NN - AFM NNN regime at U d = 2.0 and 2.5 eV. Ad- 
ditionally, the numbers are in excellent agreement with 
the TB results, the difference between J 2 and J 2 (quan- 
tified by the J 2 /J 2 ratio) as well as the different AFM 
contributions to J\ and J[ (quantified by the difference 
Jj — J 1, similar to Table [ll| are remarkably reproduced. 

The FSL-likc spin system of the AA'VO(P0 4 ) 2 com- 
pounds includes four incquivalent exchange couplings. 
In case of BaCdVO(P0 4 ) 2 , the respective NN (Ji,J[) 
and NNN (</ 2 , J 2 ) couplings nearly match and give rise 
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TABLE III: LSDA+[7 estimates of the exchange couplings in BaCdVO(P0 4 )2 and SrZnVO(P0 4 )2. U d is the Coulomb 
repulsion parameter of LSDA+f/. The columns ( J[ — Ji) and J2/J2 should be directly compared to the last column of Table [il] 



AA' 
BaCd 



SrZn 



U d (eV) 

2.0 

2.5 
3.0 

2.0 
2.5 
3.0 



Ji (K) 

-8.6 
-2.2 
1.6 

-10.0 
-2.4 
2.7 



J[ (K) 

-6.3 
-1.3 
1.6 

-5.4 
1.6 
6.3 



Ji - Ji (K) 

2.3 
0.9 
0.0 

4.6 
4.0 
3.6 



J 2 (K) 

27 
21.9 
18.1 

24.8 
20.8 
17.5 



J'2 (K) 

21.6 
17.8 
14.6 

11.3 
8.7 
7 



J2/J2 

0.80 
0.81 
0.81 

0.46 
0.42 
0.40 



to the almost regular FSL. In BaZnVO(P04)2 and 
Pb2VO(P04)2, the distortion is more pronounced. In 
case of SrZnVO(P04) 2 , the square lattice is strongly 
distorted: the two NNN couplings differ by a factor of 
three, and the two NN couplings do not match as well. 
Using J™ p ~ —8.3 K as the averaged value (Table |l| 
and assuming J[ — J x ~ 3 K (Table [TT]) , we estimate 
J[/Jl — 0.69. In the next section, we discuss the struc- 
tural origin of the spin lattice distortion. 



STRUCTURAL ORIGIN OF THE 
DISTORTION 



To study the influence of individual structural factors 
on the spin lattice distortion in the AA'VO(PC>4)2 com- 
pounds, we construct a number of model structures. The 
initial structure resembles that of ai-LiVOP04 (see Sec- 
tion . 



Ill I and includes regular vanadium and phosphorous 
polyhedra in the [VOPO4] layers (the left panel of Fig. [2]) . 
The cation-oxygen separations are 1.950 A for the basal 
oxygen atoms and 1.582 A for the axial oxygen atom 
in the VO5 square pyramids and 1.543 A for the PO4 
tetrahedra. Then, we introduce certain distortions and 
analyze their influence on the magnetic interactions. We 
use the TB approach; therefore, we focus on the AFM 
NNN couplings and trace the change of the J 2 /J 2 ratio. 
In the end of the section, we briefly comment on the NN 
couplings and the J[ vs. J\ distortion. 

The difference between the AA'VO(P04) 2 compounds 
originates from different metal cations located between 
the [VOPO4] layers. The cation size is reduced along 
the series from AA' = BaCd to AA' = SrZn, and this 
trend correlates to the reduction of J' 2 / J2 from ~ 0.8 in 
BaCdVO(P0 4 ) 2 to ~ 0.4 in SrZnVO(P0 4 )2- The coinci- 
dence of the two trends gives a hint that the cation size 
should be the origin of the distortion. Still, the influence 
of the cation size on the structure of the [VOPO4] layers 
is quite complex. 

Smaller metal cations tend to have shorter metal- 



oxygen distances and lower coordination numbers. 
Thus, barium is surrounded by nine oxygen atoms in 
BaZnVO(P0 4 ) 2 and BaCdVO(P0 4 ) 2 , while strontium 
has only eight neighboring oxygens in SrZnVO(P04)2- 
The same holds for cadmium (coordination number of 6) 
vs. zinc (coordination number of 4). Larger Ba and 
Cd cations are compatible with the nearly flat [VOPO4] 
layers (see the left panel of Fig. |6j. To provide proper 
oxygen coordination for smaller cations, the layers have 
to buckle. The Sr and Zn cations occupy the positions 
at the points of the downward and upward curvature, re- 
spectively (see the right bottom panel of Fig. [6]). The 
layer buckling can be quantified via the angle ip, see Ta- 
ble |IV| and the middle panel of Fig. [2] 

The structural changes are not confined to the layer 
buckling. In particular, the unit cell parameters change 
in a peculiar manner. As the cation size is decreased, the 
a and b parameters (in-layer spacings) are increased (see 
Table IV), while the c parameter (interlayer spacing) is 



decreased to provide the overall reduction of the unit cell 
volume expected for the substitution by a smaller cation. 
The increase of the in-layer dimensions can also be un- 
derstood via the change of the cation coordination num- 
bers. Barium cations coordinate oxygen atoms from four 
surrounding PO4 tetrahedra in the [AA'PC^] interlayer 
block (upper left panel of Fig. [6]) . Strontium cations coor- 
dinate three tetrahedra only (upper left panel of Fig. |6| , 
while the fourth tetrahedron is "pushed away", thus ex- 
panding the unit cell along the a and b directions. 

The picture presented in the last two paragraphs and 
visualized in Fig. [6] does not reflect all the structural 
changes including, e.g., slight shifts of the metal cations 
and tiltings of the PO4 tetrahedra. However, this pic- 
ture grasps the essential changes that bear influence on 
the magnetic [VOPO4] layers and on the exchange cou- 
plings. There are two main structural changes in the 
magnetic layers: (i) the buckling and (ii) the stretching 
in the ab plane. The latter is performed via the dis- 
tortion of the VO5 pyramids, while the PO4 tetrahedra 
remain rigid. The specific distortion of the square pyra- 
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FIG. 6: (Color online) Crystal structures of BaCdVO(P04)2 (left panels) and SrZnVO(P04)2 (right panels): the upper panels 
show the interlayer [AA'POzi] blocks, while the bottom panels present the buckling of the [VOPO4] layers. The change of Ba 
and Cd for Sr and Zn leads to the reduction of the coordination numbers and the resulting reorganization of the structure: the 
layer buckling and the stretching in the ab plane (see text for details). 
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FIG. 7: (Color online) Distortion of the [VOPO4] layer as 
implemented in the model structures. Arrows show displace- 
ments of the vanadium atoms away from their ideal positions 
in the centers of the VO5 square pyramids. The displacements 
yield two different V-0 distances (d and d!) and two inequiv- 
alent NNN couplings (J2 and J2). Small spheres with solid 
and hatched filling denote different types of oxygen atoms 
with shorter (white bonds) and longer (shaded bonds) V-O 
distances, respectively. Dashed and dotted lines indicate the 
interactions J2 and J 2- 



mids is shown in Fig. |7j vanadium atoms are shifted 
away from the pyramid center and yield two longer and 
two shorter V-0 bonds. Two types of the V-0 bonds are 
easily distinguished in Table[lV} The V-0(2) and V-0(8) 
distances remain nearly constant (1.95 — 2.00 A) in the 
whole AA'VO(P04)2 series, while the two other distances 



[V-0(6) and V-0(9)] are expanded from about 2.00 A in 
BaCdVO(P0 4 ) 2 up to 2.11 A in SrZnVO(P0 4 ) 2 - Both 
the buckling and the stretching of the magnetic layers 
are most pronounced in SrZnVO(P04)2- Below, we con- 
struct proper model structures and separately analyze 
the influence of these effects on the exchange couplings. 

To reproduce the layer buckling (i), we keep the VO5 
and PO4 polyhedra rigid and simply change the V-O-P 
angles to achieve the necessary buckling angle (p. To 
reproduce the layer stretching (ii), vanadium atoms are 
shifted away from the centers of the pyramids (see Fig. [7]) , 
and the unit cell is properly expanded in the ab plane. 
Then, two V-0 distances remain constant (d — 1.95 A), 
while two other distances (d f ) are increased. The expan- 
sion in quantified by the value of Ad — d' — d. 

In Fig. [8] we plot the distortion of the NNN couplings 
{J 2 /J 2) as found for the real compounds and for the 
model structures. We find that the buckling of the lay- 
ers is unable to account for the spin lattice distortion 
observed in most of the AA'VO(P04)2 compounds (see 
open circles in Fig. 111). On the other hand, the shifts of 
the vanadium atoms and the resulting layer stretching 
perfectly reproduce the distortion, even without consid- 
ering the layer buckling (see filled circles in Fig. [8j . Thus, 
we conclude that the distortion of the NNN couplings in 
AA'VO(P04)2 originates from the distortion of the VO5 
pyramids. We can also unambiguously assign the weaker 
interaction J' 2 to the longer V-0(6) and V-0(9) sepa- 
rations consistent with the LSDA+C/ results for the real 
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TABLE IV: Lattice parameters (a,b,c) and some other geometrical characteristics of the AA'VO(P04)2 compounds. The 
V-O distances are given in units of A, and the oxygen positions are numbered according to the structural data in Refs. [37| and 
1381 The angle ip is a measure for the buckling of the [VOPO4] layers as shown in Fig. p] 



AA' 


a (A) 


6(A) 


c(A) 




V-0(2) 


V-O (8) 


V-O (6) 


V-O (9) 


Ref. 


BaCd 


8.838 


8.915 


19.374 


172 


1.977 


1.975 


2.011 


1.992 


EH 


Pb 2 a 


9.016 


8.747 


9.863 6 


155 


1.954 


1.975 


2.024 


2.000 


\m 


BaZn 


8.814 


9.039 


18.538 


160 


1.956 


1.974 


2.045 


1.993 


M 


SrZn 


9.066 


9.012 


17.513 


150 


1.971 


1.999 


2.110 


2.039 





"Non-standard setting used, see Sec. [TT1 

'The unit cell of Pb2VO(P04)2 includes one magnetic layer in contrast to the other AA'VO(P04)2 compounds with two [VOPO4] layers 
in the unit cell. Yet the c parameter is not the true (shortest) interlayer distance due to the monoclinic symmetry of the structure. 
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FIG. 8: (Color online) Spin lattice distortion (J2/J2) for 
the model structures with the layer buckling (empty cir- 
cles) and the layer stretching (filled circles) and for the real 
AA'VO(P04)2 compounds (filled diamonds). The layer buck- 
ling is quantified by the buckling angle tp (Fig. . The layer 
stretching is imposed by shiftingvanadium atoms and quan- 
tified by Ad = d! — d (see Fig. [7b. For the real compounds, 
d = (dv-o(a) +d V -o(8))/2 and d = (d V -o(e) + <*v-o(9))/2 (see 
Table pV]. 



compounds. 

Unfortunately, the trends for the NN couplings are less 
clear. According to the discussion in Sec. |IV[ the differ- 
ence between Ji and J[ originates from different AFM 
contributions to these couplings. The initial model struc- 
ture does not show any sizable NN hoppings (both t\ 
and t'i are below 4 meV), and the shifts of the vana- 
dium atoms within the flat layer do not change these 
hoppings. The layer buckling enlarges t\ and t\ up to 
10 - 15 meV, i.e., the TB results for BaZnVO(P0 4 ) 2 are 
reproduced (see Table |n|. However, the model structures 
do not show the anisotropy of the NN couplings, as found 
in SrZnVO(P0 4 ) 2 and Pb 2 VO(P0 4 ) 2 . The anisotropy is 
likely caused by more subtle changes that are not in- 
cluded in our model structures. Nevertheless, the flat 



[VOPO4] layer [as found in BaCdVO(P0 4 ) 2 ] does not 
show any considerable AFM contributions to the NN cou- 
plings, hence the anisotropy of the NN couplings in the 
flat-layer compounds should be small. 



VI. EXTENDED FSL MODEL 

In this section, we discuss the properties of the ex- 
tended FSL model. We address thermodynamic prop- 
erties, magnetic susceptibility and specific heat, because 
these quantities are measured experimentally and com- 
monly used for the evaluation of the exchange couplings 
in the FSL compounds. The extended model includes 
four independent parameters, but we are mainly inter- 
ested in the role of the distortion, i.e., the difference be- 
tween J 1 and J[ or J 2 and J 2 . Therefore, we fix the av- 
eraged NN and NNN couplings (Ji and J 2 , respectively) 
and the effective frustration ratio a = 3%jJ\. Then we 
vary either J\ and J{ or J 2 and J 2 . We performed the 
simulations for two representative values, a = —2 and 
a = —1, to study the frustration regime relevant for 
Pb 2 VO(P0 4 ) 2 and BaZnVO(P0 4 ) 2 or SrZnVO(P0 4 ) 2 
and BaCdVO(P0 4 ) 2 , respectively (see Table [[]). In the 
following, we present the results obtained at a = —2 
(Fig. [9]). The simulations for a = — 1 reveal a very simi- 
lar behavior, thus we do not discuss them in detail. 

The results for the regular FSL at a = —2 match that 
of Ref. [HI (for the comparison, one should use the frus- 
tration angle tyf/it — 0.65 with ipf — tan" 1 ^). The 
distortions of the NN and NNN bonds have different ef- 
fects on the thermodynamic properties. The distortion 
of the NNN bonds ( J 2 vs. J 2 , left panel of Fig. ^ leads 
to a slight shift of the susceptibility maximum to lower 
temperatures, while the absolute value at the maximum 
is increased. The maximum of the specific heat is also 
shifted to lower temperatures, but its height is reduced. 
The distortion of the NN bonds (J[ vs. J\, right panel of 
Fig. [9| leads to opposite and more pronounced changes. 
The maxima are shifted to higher temperatures, the sus- 
ceptibility maximum is reduced, while the specific heat 
maximum is increased. 
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FIG. 9: (Color online) Full diagonalization results for the distorted FSL model: magnetic susceptibility (primary figures) and 
specific heat (insets). The simulations are performed at fixed averaged couplings Ji = ( Ji + J[)/2 and J2 = (J2 + J^)/2 and the 
fixed frustration ratio J2/J1 = — 2. In the left and right panels, J2 and J'2 or Ji and J[ are varied, respectively. The arrows show 
the changes upon increasing the distortion. The thermodynamic energy scale J c is defined as J c = \J ( J\ + J{ 2 + J2 + Jj 2 )/2. 



To understand these results, one should recall that the 
magnetic susceptibility and the specific heat maxima are 
caused by correlated spin excitations. These excitations 
compete with quantum fluctuations caused by the low 
dimensionality and the magnetic frustration. Thus, the 
shift of the maxima to lower/higher temperatures implies 
the enhancement/reduction of the quantum fluctuations. 
Then the effects for the susceptibility and the specific 
heat are opposite, since spin correlations contribute to 
the specific heat and increase its value, but lead to AFM 
ordering, hence reducing the susceptibility. The trends 
presented in Fig.[9]suggest that the distortion of the NNN 
bonds enhances quantum fluctuations (see the explana- 
tion below), while the distortion of the NN bonds reduces 
the fluctuations, and the latter effect is more pronounced. 

The effect of the NN bonds distortion is consistent 
with the previous the oretic al results for the spatially 
anisotropic FSL modelP" ^ 12 * 14 ! The narrowing and clos- 
ing of the critical (spin-liquid) region corresponds to the 
reduction of the quantum fluctuations, as observed in the 
thermodynamic data. One can get further insight into 
this effect by considering energies of the ordered struc- 
tures within the classical Heisenberg model. At a = —2, 
the competing ground states are the FM and columnar 
AFM ordering (see Fig. [TJ. The columnar AFM state is 
favored by AFM NNN interactions and by the FM in- 
teraction along the directions of columns (say, along the 
b axis, i.e., the respective interaction is J\). Yet, the 
FM interaction J[ along the a axis is unfavorable for the 
columnar ordering. As the absolute value of J[ is reduced 
and that of J\ is increased, the columnar AFM state is 
stabilized, and frustration is released. Applying similar 
considerations to the NNN bonds distortion, we find that 
the distortion does not change the energies of the FM 
and columnar AFM states, hence the magnitude of the 
frustration should remain unchanged. This conclusion 



is consistent with the relatively weak effect of the NNN 
bonds distortion. Still, the certain enhancement of the 
frustration is clearly visible in the thermodynamic data 
and likely related to quantum effects. We can speculate 
that the reduction of the J' 2 / J2 ratio leads to the forma- 
tion of spin chains within the 2D lattice. The chains are 
formed by the J 2 bonds and run along the b axis (see the 
right panel of Fig.|2|. Then, this one-dimensional feature 
of the spin system should enhance quantum fluctuations 
due to the effectively reduced dimensionality. 

Now, we turn to the case of moderate distortion 
relevant for the AA'VO(PC>4)2 compounds. Accord- 
ing to Sec. |IV| the strongest distortion is found in 
SrZnVO(P0 4 ) 2 with J[/J x ~ 0.7 and J' 2 /J 2 ~ 0.4. 
The respective susceptibility curves in Fig. [9] nearly 
match that for the regular FSL. Thus, the fitting of 
the experimental susceptibility data should yield aver- 
aged exchange couplings of the distorted FSL, Ji and 
J 2 . This conclusion provides a reliable basis for the in- 
terpretation of the experimental values listed in Table [T] 
The changes in the specific heat curves are also minor, 
hence the experimental specific heat will be described 
by the HTSE. Thus, the thermodynamic properties of 
the AA'VO(P04)2 compounds should fit the regular FSL 
model with averaged exchange para meters a s an excellent 
approximation, and this is the case!2412S12QI 



VII. DISCUSSION AND SUMMARY 

In this study, we performed a detailed microscopic 
investigation of the distorted FSL spin systems in the 
AA'VO(P04)2 compounds. We estimated the magnitude 
of the distortion, found the structural origin of the dis- 
tortion and analyzed the thermodynamic properties of 
the extended FSL model. Below, we consider the con- 
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sequences of these findings in several aspects: (i) the 
physics of the AA'VO(P04)2 phosphates; (ii) layered 
vanadium compounds as a playground for the search of 
new FSL materials; and (iii) lattice distortion in frus- 
trated spin systems. 

The basic experimental results for the AA'VO(PC>4)2 
compounds can be int erpreted within the fra mewo rk of 
the regular FSL model.^ 25 * 30 1 According to Sec.[Vl] ther- 
modynamic measurements lead to relevant, averaged ex- 
change couplings J° xp and J%* p . These values place all 
the AA'VO(PC>4)2 compounds to the columnar AFM re- 
gion of the FSL phase diagram (Fig. [I]). Indeed, neutron 
scattering results for Pb 2 VO(P0 4 ) 2 and SrZnVO(P0 4 ) 2 
(Refs. [26] and I27p confirm the columnar ordering. Yet, 
the experimental situation is not fully clear, since nuclear 
magnetic resonance (NMR) and muon spin relaxation 
(/iSR) studies suggest a broad distribution of local mag- 
netic fields in the ordered phase, hence pointing to a more 
complex ground state (e.g., incommensurate ordering), at 
least for some of the AA'VO(PC>4)2 compounds.^ One 
can suggest that the ground state of these materials is 
influenced by the spin lattice distortion. Although de- 
tailed investigation of the ground state of the extended 
FSL model lies beyond the scope of the present study, 
the introduction to the extended model and the evalua- 
tion of the model parameters is a first step toward un- 
derstanding the ground state of AA'VO(P04)2- Further 
theoretical and experimental (neutron scattering, NMR, 
/iSR) studies on the ground state properties are highly 
desirable and should be stimulated by our work. 

Apart from the thorough studies of the ground state, 
one can suggest a more simple way for the exper- 
imental observation of the spin lattice distortion 
AA'VO(P0 4 ) 2 - According to Sec 
of the nearest-neighbor bonds (Ji 



VI 



vs. 



m 

the distortion 
J[) stabilizes 

the columnar AFM state with respect to the FM state. 
Within the classical model, the energy difference be- 
tween the FM and columnar AFM states corresponds 
to the saturation field. Therefore, the J\ ^ J[ scenario 
should have an effect on the saturation field. In case of 
BaCdVO(PC>4)2 with the nearly regular FSL, the satura- 
tion field is in excellent agreement with the averaged cou- 
plings J^ xp and (Ref . ISD]) . However, the saturation 
field of SrZnVO(P0 4 ) 2 should be different from the field 
estimated using J° xp and J^- High- field magnetization 
studies of the AA'VO(P0 4 )2 compounds will challenge 
this proposition and enable the quantitative analysis of 
the spin lattice distortion in AA'VO(P04)2- 

The next important issue is the capability of layered 
vanadium compounds to reveal new strongly frustrated 
FSL materials. The basic structural clement of these 
compounds is the [VOXO4] layer shown in the left panel 
of Fig. [2] In this layer, non-magnetic X0 4 tetrahedra 
mediate either FM or weak AFM NN and AFM NNN in- 
teractions. Then there are two ways to reach the strongly 
frustrated regime of a ~ —0.5: one should cither increase 
the absolute value of Ji [as observed in SrZnVO(PC>4) 2 ] 
or decrease J2 [as observed in BaCdVO(P04)2]- Accord- 



ing to the results of our study (Sec. [V}, the NNN inter- 
actions (J2) are sensitive to the V-0 distances, since the 
magnitude of any superexchange interaction depends on 
the orbital overlap. Thus, it is possible to reduce the 
J 2 value by increasing the V-0 separations in the VO5 
square pyramids. The factors influencing on the value of 
Ji are less clear. In the AA'VO(P04)2 compounds, the 
Ji values correlate with the distortion of the NNN bonds. 
One may suggest the layer buckling or the distortion of 
the VO5 pyramids as possible reasons for the increase of 
J° xp from BaCdVO(P0 4 ) 2 to SrZnVO(P0 4 ) 2 . However, 
this conclusion is rather empirical, and other structural 
factors may be relevant as well. 

The proper FSL material should combine the strong 
frustration with the lack of the distortion or, at least, 
with a relatively weak distortion of the spin lattice. It is 
also desirable to find a family of isostructural FSL com- 
pounds. Then, the replacement of the metal cations may 
facilitate the tuning of the system toward the strongly 
frustrated regime. Our results provide a clear recipe 
for the search of undistorted FSL materials. To get a 
regular FSL, one should keep the magnetic layer flat 
and avoid the distortion of the VO5 square pyramids. 
Clearly, it is quite difficult to fulfill these criteria within 
the AA'VO(P04)2 series. Different metal cations require 
different coordination numbers (see Fig. [6]); therefore, 
most of the respective compounds reveal a distorted FSL. 
In BaCdVO(P04)2, the cation sizes are optimal to yield 
nearly flat and regular magnetic layers and to result in 
the weak distortion of the FSL. Though tuning is pos- 
sible (see Table |l| , most of the metal cations lead to a 
layer distortion and to a spin lattice distortion as well. 

To avoid the spin lattice distortion within a com- 
pound family, one can try to reduce the number of 
metal cations and to look for another filler of the in- 
terlayer space. For example, one can consider the 
A ( VOPO4) 2 ; 4H 2 O compounds (A = Ca, Cd, Sr, Pb, Ba, 
and Mg/ZnP ^" ' 01 '^ ' 63 ' 64 ' ^ that reveal layered struc- 
tures with [VOPO4] layers separated by metal cations 
and water molecules. The resulting symmetry is or- 
thorhombic, monoclinic, or even triclinic, but the respec- 
tive spin lattice distortion should be weak. The layers are 
nearly flat for all the six compounds reported, and the 
distortion of the VO5 pyramids is also negligible. Yet, the 
change of the metal cation bears influence on the mag- 
netic interact ions, as indi cated by different Curie- Weiss 
temperatures (SSE2ME1M1 The magnitude of the frustra- 
tion in the A(VOP0 4 ) 2 • LH2O compounds remains un- 
known. All the reports available consider NN interac- 
tions only. Clearly, this scenario is oversimplified, and 
one has to apply the FSL (rather than a simple square 
lattice) model while analyzing magnetic properties of 
A(VOP04)2 • 4H 2 0. Further studies of these materials 
should be very promising. 

The last, but not least, point deals with the influence 
of the distortion on the properties of frustrated spin sys- 
tems. In this study, we exemplified this problem by con- 
sidering the distortion of the FSL. We found that differ- 
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ent types of distortion had different effects on the magni- 
tude of the frustration. The distortion of the NN bonds 
stabilizes the columnar AFM ordering and releases the 
frustration. Yet, the distortion of the NNN bonds keeps 
the strong frustration and even enhances the quantum 
fluctuations. These results suggest that both regular 
and distorted frustrated materials should be considered 
as proper realizations of the frustrated spin models. The 
case of the AA'VO(PC>4)2 compounds shows how the dis- 
torted FSL materials can be successfully treated within 
the regular FSL model. In these compounds, the ther- 
modynamic properties at available temperatures nearly 
match those of the regular model. However, the ground 
state properties may be different. We believe that the 
study of frustrated materials with spin lattice distortion 
will improve our understanding of the frustrated spin sys- 
tems and facilitate the observation of interesting phenom- 
ena suggested by theory. 

In conclusion, we have studied the distorted frus- 
trated square lattice in layered vanadium phosphates 
AA'VO(P0 4 ) 2 (AA' = Pb 2 , SrZn, BaZn, and BaCd). 
In these compounds, both nearest-neighbor and next- 
nearest-neighbor bonds of the square lattice are dis- 
torted, hence an extended spin model with four inequiv- 
alent exchange couplings should be considered. Our es- 
timates of the individual model parameters suggest the 
least pronounced distortion in BaCdVO(PC>4)2 and the 
most pronounced distortion in SrZnVC^PO^. The dif- 
ference between the nearest-neighbor and next-nearest- 
neighbor interactions in these compounds originates from 



peculiar structural changes upon substituting interlayer 
metal cations A and A'. The buckling of the magnetic 
layers may change the interactions along the side of the 
square, while the distortion of the vanadium polyhedra 
and the resulting stretching of the magnetic layer lead 
to the difference between the diagonal interactions of the 
distorted square lattice. The distortion of the square 
lattice in AA'VO(PC>4)2 is moderate from the point of 
view of the thermodynamic properties. The tempera- 
ture dependences of the magnetic susceptibility and the 
specific heat resemble those of the regular square lattice, 
hence previous experimental studies of AA'VO(PC>4)2 re- 
ported averaged couplings. These couplings can be used 
to place the compounds on the phase diagram of the reg- 
ular model. In contrast to the moderate influence on 
the thermodynamic properties, the spin lattice distortion 
may have a larger effect on the ground state. Further ex- 
perimental and theoretical studies of this problem are 
highly desirable. 
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